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Abstract 

Theoretical reaction rates in the temperature range 0.01 x 10 9 < T(K) < 
10.0 x 10 9 are calculated in the statistical model (Hauser-Feshbach formal- 
ism) for targets with 10 < Z < 83 (Ne to Bi) and for a mass range reaching 
the neutron and proton driplines. Reactions considered are (n,7), (n,p), 
(n,a), (p,7), (p,a), (a, 7), and their inverse reactions. Reaction rates as 
a function of temperature for thermally populated targets are given by 
analytic seven parameter fits. To facilitate comparison with experimental 
rates, the stellar enhancement factors are also tabulated. Two complete 
sets of rates have been calculated, one of which includes a phenomeno- 
logical treatment of shell quenching for neutron-rich nuclei. These ex- 
tensive datasets are provided on-line as electronic files, while a selected 
subset from one calculation is given as printed tables. A summary of the 
theoretical inputs and advice on the use of the provided tabulations is 
included. 
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1 Introduction 

Nuclear reaction rates are an essential ingredient for all investigations of nu- 
cleosynthetic or energy generating processes in astrophysics. Highly unstable 
nuclei are produced in such processes which again can become targets for sub- 
sequent reactions. Cross sections and astrophysical reaction rates for a large 
number of nuclei are required to perform complete network calculations which 
take into account all possible reaction links and do not postulate a priori sim- 
plifications. Despite concerted experimental efforts, most of the involved nuclei 
are currently not accessible in the laboratory and therefore theoretical models 
have to be invoked in order to predict reaction rates. 

In predictions of cross sections and reaction rates for astrophysical appli- 
cations, slightly different points are emphasized than in pure nuclear physics 
investigations. Firstly, one is confined to the very low-energy region, from ther- 
mal energies up to a few MeV. Secondly, since most of the ingredients for the 
calculations are experimentally undetermined, one has to develop reliable phe- 
nomenological models to predict these properties with an acceptable accuracy 
across the nuclear chart. This task is made even harder by the lack of infor- 
mation on specific properties, such as optical potentials for a particles at the 
astrophysically relevant energies even for stable nuclei. Therefore, one has to be 
satisfied with a somewhat more limited accuracy, as compared to usual nuclear 
physics standards. The accuracy of the rates is estimated to be within a factor 
of 1.5 — 2, with an even better average deviation, e.g. of 1.4 for neutron capture. 
Considering the substantially larger uncertainties in many astrophysical scenar- 
ios, this seems to be acceptable. Thus, the real challenge is not the application 
of well-established models, but rather to provide all the necessary ingredients 
in as reliable a way as possible, also for nuclei where no such information is 
available. Many efforts have been directed at addressing those problems and 
the current status of the investigations make it worthwhile to publish a full 
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set of theoretical rates, intended to supersede early reaction rate tabulations 

DBiiii. 

For the majority of nuclear reactions in astrophysics, the statistical model 
(Wolfenstein-Hauser-Feshbach approach) [Q, ^ can be applied. This is appro- 
priate provided the level density in the contributing energy window around the 
peak of the projectile energy distribution is sufficiently high to justify a statis- 
tical treatment. The critical level density is usually estimated between 5 and 
10 MeV" 1 ||. Furthermore, the compound nucleus picture will only dominate 
when the energy of the incident particle is low enough (< 20 MeV). While the 
latter point is practically always satisfied in astrophysical environments, the 
level density may fall below the critical value in certain nuclei lighter than Fe, 
at shell closures, and for very neutron-rich or proton-rich isotopes near the drip 
lines with correspondingly low separation energies. In these cases, single reso- 
nances or direct capture contributions will become significant and have to be 
treated individually. In this tabulation, we are not concerned with such effects 
but rather give a full set of rates calculated in the stat istical model. However, 
the limits of its applicability will be discussed (Sec. 2.4). 

In the following section, we concisely summarize the theoretical background 
for easy reference as well as the nuclear properties used as input in the calcu- 
lations. This is followed by a section defining the reaction rates, explaining the 
fitting procedure and giving more details on the tabulated values. The paper is 
concluded by a summary and the rate tables. 



2 The Statistical Model 
2.1 Theory 



The averaged transmission coefficients T comprise the central quantities in sta- 
tistical model calculations. They do not reflect a resonance behavior, but rather 
describe absorption via an imaginary part in the (optical) nucleon-nucleus po- 
tential [fto"[ . This leads to the well-known expression 

{Eij) ~ (2jr+D(2j J+ i) (1) 

Tf{E, J, jT, E%, J?,*?)TZ(E, J, tt, r m , tO 

x 5 (2J + 1} 

for the cross section a^ v of the reaction i^(j,6)m v from the target state z M to 
the excited state m v of the final nucleus, with a center of mass energy E^j and 
reduced mass puj. J denotes the spin, E the corresponding excitation energy, 
and 7r the parity of excited states. When these properties are used without 
subscripts they describe the compound nucleus; subscripts refer to states of the 
participating nuclei in the reaction i^(J, o)m v and superscripts indicate the spe- 
cific excited states. The total transmission coefficient Ttot = Ylu o describes 
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the transmission into all possible bound and unbound states v in all energeti- 
cally accessible exit channels o (including the entrance channel i). Experiments 
measure er lab =Ylu ^"(-^i)) summed over all excited states of the final nucleus, 
with the target in the ground state. Target states p in an astrophysical plasma 
of temperature T* are thermally populated and the astrophysical cross section 
a* is given by 



, _ E,(2Jf + l)exp(-ff,r/fcT*)E^(^) 
a[ij) ~ V(2Jf + l)exp(- J B.f/fcr*) ' U! 



k being the Boltzmann constant. The summation over v replaces T%(E, J, 7f) in 
Eq. (|l|) by the total transmission coefficient 

T (E,J,n) = J2T^E,J,n,E^,J^n v m ) (3) 

i/=0 

rE-S m , 

+ / /J T {E, J,ir, E m , J m ,ir m )p(E m , J m ,ir m )dE m . 

Here S m . is the channel separation energy, and the summation over excited 
states above the highest experimentally known state v m is changed to an inte- 
gration over the level density p. The summation over target states p in Eq. (||) 
has to be generalized accordingly. 

The important ingredients of statistical model calculations as indicated in 
Eqs. (|l|) through (§) are the particle and 7-transmission coefficients T and the 
level density of excited states p. Therefore, the reliability of such calculations 
is determined by the accuracy with which these components can be evaluated 
(often for unstable nuclei) . It is in these quantities that various statistical model 
calculations differ. The reaction rates given in this paper are calculated with 
the code NON-SMOKER Q, derived from the well-known SMOKER code §]. 
(The code MOST Jlf] is another code derived from SMOKER.) In the following 
we want to briefly outline the methods utilized in that code to estimate those 
nuclear properties. The challenge is in the goal to provide them in as reliable a 
way as possible, also for unstable nuclei for which no experimental information 
is available. Thus, global descriptions are employed which minimize the overall 
error and are trusted to be reliable also far from stability. 



2.2 Transmission Coefficients 

The transition from an excited state in the compound nucleus (E 7 J, ir) to the 
state (E-\ Jf , 7rf ) in nucleus i via the emission of a particle j is given by a 
summation over all quantum mechanically allowed partial waves 

J+s j^+Jj 

T?(E,J,n,E?,J?,n!t)= £ E T ii.( E &> ( 4 ) 

l=\J-s\ s=\J%-Jj\ 
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Here the angular momentum / and the channel spin s = Jj + Jf couple to 
J = 1 + s. The transition energy in channel j is E^-=E — Sj — where Sj is 
the channel separation energy. 

The total transmission coefficients for this tabulation are then calculated by 
applying Eq. (||) and utilizing up to 19 experimentally known excited states. 
The data are taken from Jl]|, up to the first level for which the spin assignment 
was not known. Ground state spin and parities are known for many unstable 
nuclei. Far off stability, ground state spins and parities are taken from |l^], if 
experimental values are not available. 



2.2.1 Particle Transmission Coefficients 

The individual particle transmission coefficients Tj u are calculated by solving 
the Schrodinger equation with an optical potential for the particle-nucleus inter- 
action. We employ the optical potential for neutrons and protons given by |fL5| , 
based on microscopic infinite nuclear matter calculations for a given density, 
applied with a local density approximation. It includes corrections of the imag- 
inary part flU, 0. 

The optical potential for a particles from was shown to be quite accurate 
for a wide range of nuclei and is used in this work. However, it was realized 
[ ^9| |2(i[| that for heavily charged nuclei a more sophisticated potential had to be 
adopted at the comparatively low energies of astrophysical interest. Promising is 
the folding approach J2l[], with a parameterized mass- and energy-dependence of 
the real volume integral (22). Microscopic and deformation information should 
be considered in the parametrization of the imaginary potential ^3|. However, 
due to the scarcity of experimental data, the potential parameters can as yet 
only be extracted for a limited mass and energy range. The optical a+nucleus 
potential is likely to introduce the largest uncertainties in the charged particle 
rates presented here. Further experimental work is clearly necessary and most 
welcome. 

Deformed nuclei are treated by using an effective spherical potential of equal 
volume, based on averaging the deformed potential over all possible angles be- 
tween the incoming particle and the orientation of the deformed nucleus. 



2.2.2 Radiative transmission coefficients 

At least the dominant 7-transitions (El and Ml) have to be included in the 
calculation of the total photon width. The smaller, and therefore less important, 
Ml transitions are treated, as usual, in the simple single particle approach 
(T oc E 3 |2^|), as also discussed in The El transitions are calculated on the 
basis of the Lorentzian representation of the Giant Dipole Resonance (GDR). 
Within this model, the El transmission coefficient for the transition emitting a 
photon of energy E 1 in a compound nucleus is given by 

T (F] _8NZe*l + X ^i r Gii E* 

lEl(E y ) - -— — — — -j- > „ N2 , T-2 S3 ' i b > 



3 A ncMc 2 ^ 3 (^_^. )2 + r 2 G ^ 
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Here, M is the proton mass, \(— 0.2) accounts for the neutron-proton exchange 
contribution J25|] , and the summation over i includes two terms which correspond 
to the split of the GDR in statically deformed nuclei, with oscillations along 
(i = 1) and perpendicular (i = 2) to the axis of rotational symmetry. Many 
microscopic and macroscopic models have been devoted to the calculation of 
the GDR energies (Eg) and widths (Tg)- Here, the (hydrodynamic) droplet 
model approach |^6| is used for Eg, which gives an excellent fit to the GDR 
energies and can also predict the split of the resonance for deformed nuclei, when 
making use of the deformation, calculated within the droplet model. In that 
case, the two resonance energies are related to the mean value calculated by the 
relations @ E GA + 2E G , 2 = 3E G , E G a/E G ,i = 0.91 177 + 0.089. 77 is the ratio of 
the diameter along the nuclear symmetry axis to the diameter perpendicular to 
it, and is obtained from the experimentally known deformation or mass model 
predictions. For the width Tg of the GDR the description of is used, which 
applies to spherical and deformed nuclei and can be described as a superposition 
of a macroscopic width due to the viscosity of the nuclear fluid and a coupling 
to quadrupole surface vibrations of the nucleus (see also ||). 

Direct application of Eq. (||) would overestimate the radiation width by 
about 30% (see e.g. [^9[ Q). This is due to the fact that, for low energy 
7-transitions, the Lorentz curve is suppressed and the GDR width increases 
with excitation energy (e.g. (3l], [32)). To account for these deficiencies, various 
treatments of an energy-dependent width have been suggested. We use the 
form H 

t g {e 1 ) = t g ^l. (6) 

Another effect has to be taken into account for certain a-capture reactions. 
Because of isospin selection rules, 7-transitions between levels with isospin 1 = 
are forbidden. This leads to a suppression of the cross section for (a,j) reactions 
on self-conjugate (N = Z) targets, due to isospin conservation. A suppression 
could also be found for capture reactions leading into self-conjugate nuclei, al- 
though somewhat less pronounced because 1=1 states can be populated ac- 
cording to the isospin coupling coefficients. The suppression is usually treated 
as a suppression of the 7-width. In previous rate tabulations it was either ne- 
glected H or accounted for in a phenomenological way by dividing the 7-width 
by quite arbitrary factors [§, §|. In the code NON-SMOKER the appropriate 
7-widths are obtained by explicitly accounting for isospin mixing and suppres- 
sion of the appropriate 7-transitions |fj"l|| . A detailed account of the procedure 
can be found in f^jl - 

2.2.3 Width fluctuation corrections 

In addition to the ingredients required for Eq. (|l]), like the transmission co- 
efficients for particles and photons, width fluctuation corrections W(j, o, J, n) 
have to be employed. They define the correlation factors with which all partial 
channels for an incoming particle j and outgoing particle o, passing through the 
excited state (E, J,ir), have to be multiplied. This takes into account that the 
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decay of the state is not fully statistical, but some memory of the way of forma- 
tion is retained and influences the available decay choices. The major effect is 
elastic scattering, the incoming particle can be immediately re-emitted before 
the nucleus equilibrates. Once the particle is absorbed and not re-emitted in the 
very first (pre-compound) step, the equilibration is very likely. This corresponds 
to enhancing the elastic channel by a factor Wj . In order to conserve the total 
cross section, the individual transmission coefficients in the outgoing channels 
have to be renormalized to Tj. The total cross section is proportional to Tj 
and, when summing over the elastic channel (WjTj) and all outgoing channels 
- Tj), one obtains the condition T^T'^W^/T^) + Tj(T/ oi - Tj)/T/ ot 
We can (almost) solve for Tj 

T- 

T' = 3 - . (7) 

> l + Tj(W/,-l)/T/ ot {> 

This requires an iterative solution for T' (starting in the first iteration with Tj 
and T to t), which converges rapidly. The enhancement factor Wj has to be known 
in order to apply Eq. (Q) . A general expression in closed form was derived j34| , 
but is computationally expensive to use. A fit to results from Monte Carlo 
calculations gave |3^] 

^ = 1 + T7^ • («) 

For a general discussion of approximation methods see J36| , [JTJ. Eqs. (^) and 
(^|) redefine the transmission coefficients of Eq. (|l|) in such a manner that the 
total width is redistributed by enhancing the elastic channel and weak channels 
over the dominant one. Cross sections near threshold energies of new channel 
openings, where very different channel strengths exist, can only be described 
correctly when taking width fluctuation corrections into account. The width 
fluctuation corrections of [[35) are only an approximation to the correct treat- 
ment. However, it was shown that they are quite adequate |j38|. 



2.3 Level Densities 

Until recently, the nuclear level density has given rise to the largest uncertainties 
in the description of nuclear reactions in the statistical model [§, [|, ||. Imple- 
mented in the NON-SMOKER code is a recently improved treatment |J. It 
is based on a shifted Fermi-gas formalism |39| with an energy-dependent level 
density parameter a together with microscopic corrections from nuclear mass 
models. This leads to improved fits to known level densities in a wide range of 
masses More sophisticated Monte Carlo shell model calculations as well 
as combinatorial approaches (see e.g. 0]), have shown excellent agreement with 
this phenomenological approach and justified the application of the Fermi-gas 
description at and above the neutron separation energy. 

An in-depth description of the model and its application to astrophysical 
problems can be found in ||. Here, we only want to briefly summarize the inputs 
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used for calculating the rates presented in this tabulation. It should be noted 
that we applied our description throughout the nuclear chart, without relying on 
experimental level density parameters in specific cases as has been done before 
[|[ ||, This may lead to locally slightly larger deviations from experiment 
but it improves the reliability when extrapolating to unknown isotopes. 

The microscopic correction and the pairing corrections comprise crucial in- 
puts for the level density formalism used here (see J^] for details). They can 
be extracted from mass models. There is a choice of several mass models in 
the NON-SMOKER code. The Finite Range Droplet Model (FRDM) @ and 
an extended Thomas-Fermi approach with Strutinski Integral (ETFSI-Q) |43| 
have been chosen for the reaction rate calculations in this work. It has to be 
emphasized that experimental mass values JZ^] were included where possible. 
This is straightforward for the separation energies which were calculated from 
the mass differences; it was ensured that either only experimental or theoretical 
values were used in the differences, thus avoiding unphysical breaks at transition 
points from experiment to theory. The microscopic corrections were obtained by 
subtracting the well-defined spherical macroscopic (droplet) term of the FRDM 
from the total mass energy derived from experiment, from the FRDM or from 
ETFSI-Q, respectively (cf. Eq. (17) in ||). The validity of the resulting rates is 
discussed in Sec. Rates based on other mass models can be obtained from 
the authors on request or on-line (see Sec. ^). 

The shifted Fermi-gas approach diverges for U = E — 8 = (i.e. E = 8, 
if 8 is a positive backshift, with E being the excitation energy and 8 being an 
energy shift due to pairing corrections) . In order to obtain the correct behavior 
at very low excitation energies, the Fermi- gas description can be combined with 
the constant temperature formula [Q and references therein) 

P (u) oc . (9) 

nucl 

The two formulations are matched by a tangential fit determining the nuclear 
temperature T nuc i. 



2.4 Applicability of the Statistical Model 

The statistical model can be applied provided that the use of averaged transmis- 
sion coefficients (Eq. (^)) is permitted. This will be the case for high level den- 
sities with completely overlapping resonances, typical for the compound nucleus 
reaction mechanism. For light nuclei, decreasing particle separation energies or 
at shell closures, level densities will eventually become too low for the appli- 
cation of the statistical model at astrophysical temperatures. In those cases, 
single resonances and contributions from the direct reaction mechanism have to 
be t aken into account jlH]]. Based on the level density description outlined in 



Sec. [2.3| , a quantitative criterion for the applicability was derived recently 
In the present work we give tables of all reaction rates regardless of applicability 
but specify the allowed temperature range in the tables. The estimate is quite 
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conservative and thus the rates can still be accurate slightly below the given 
lower limits of the temperature. 



3 Astrophysical Reaction Rates 
3.1 Definition 

The nuclear reaction rate per particle pair at a given stellar temperature T* is 
determined by folding the reaction cross section cr*(E) from Eq. (||) with the 
Maxwell-Boltzmann velocity distribution of the projectiles ]46| 

1/2 00 

«"•■>-«">•-(£) ^h E ) E M:<P) iE ■ < 10 > 

It has to be emphasized that only the use of the stellar cross section a* (Eq. (g)) 
yields a reaction rate with the desired behavior that the inverse reaction can 
be calculated by using detailed balance. Therefore, laboratory rates - which 
only measure cr lab = ~}2 v <J au , i.e. the cross section with the target being in 
the ground state - should always be measured in the direction that is least 
affected by excited target states. This is usually the reaction with positive Q- 
value (exoergic reaction). For astrophysical applications, such rates have to be 
corrected for the stellar enhancement effect due to the thermal excitation of the 
target [Q. The stellar enhancement factors (SEF) /* are defined by 

r = £s ■ ( n ) 

The values of /*for a range of temperatures for nuclei close to stability are given 
in the tables. Stellar enhancement factors for neutron capture reactions and a 
discussion of the involved uncertainties can also be found in a recent compilation 
of neutron cross sections for the s process p7| . 



3.2 Partition functions and reverse rates 

The temperature-dependent partition function G(T*) normalized to the ground 
state target spin Jf is defined as J|^| 

Urn 

(2J° + 1)G(T*) = ^(2Jf + l)e' E "' kT ' (12) 

+ / E (2J /i + l)e- e/feT *pM / ^)* , 

with p being the level density and \i m the last included experimentally known 
state. The included experimental levels were the same as for the calculation of 
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the transmission coefficients (Eq. (||)). For the temperature range considered 
here, the maximum energy J EP iax above which there are no more significant 
contributions to the partition function is of the order of 20 — 30 MeV. With 
that definition, the stellar reaction rate (<iiv)* for a reaction with particles in 
all channels is related to the rate of the reverse reaction (a m v)* by 

N A (a y)' ( A ^ (2J - + 1)(2 ^ ± 1} e~^N A (a v)* 

{ m ) ~ \A a A m ) (2J + l)(2J m + 1) G m (T*) E NA {aiV) ■ 

(13) 

where Na is Avogadro's number, J and A are spins and masses A (in atomic 
mass units u) of the particles involved in the reaction i(j,o)m, Q is the re- 
action Q-v&\ue. To calculate photodisintegration rates from capture rates the 
appropriate relation is 

_ (AAA 3/2 (2J i + l)(2J J + l) Gj(T*) 3/2 Q/kT * , 

X ^Va^) (2J m + l) G m {T*) iT) Fe NaM ■ 

(14) 

For A 7 in s 1 and using the usual practical units, i.e. temperatures Tg = T*/10 9 
K and Na ((Tv) in cm 3 s _1 mole -1 , one obtains 

(T*) 3/2 F = ( = T 9 3/2 9.8685 x 10 9 mole cm" 3 . (15) 
V ; \2nti 2 J N A 9 K J 

The numerical factor F as well as the spin and mass factors are already ac- 
counted for in the parameter cig cv which is tabulated. See Sec. 3.3.2) for further 
details. 



In the tabulated rates, the thermal effects are already considered. Statisti- 
cal model calculations not including full photon cascades may also be prone to 
some error arising from the decay of unbound particle states in reactions with 
negative Q- va l ues (endoergic reactions). Furthermore, due to the exponential 
dependence of the inverse rate on the Q-value (Eqs. ( p~3| ) and (p^|)), inaccura- 
cies in the rate would be strongly enhanced when computing the exoergic rate 
from the endoergic one. In order to minimize the error, reactions are always 
calculated in the exoergic direction (with the exceptions of capture reactions, 
where photodisintegration is always treated as the reverse reaction regardless of 
Q-value) and detailed balance is applied to obtain the values for the endoergic 
reaction. This treatment has the additional advantage that it ensures consis- 
tent values for forward and reverse reactions, which is essential for application 
in astrophysical nuclear reaction networks. To calculate the actual (endoer- 
gic) reaction rate, those fits have to be multiplied by the ratio of the partition 
functions of the final nucleus and the target Gi/G m at the appropriate temper- 



ature (see Sec. 3.3.2 and Sec. 3.3. (.). For that purpose, partition functions are 



tabulated separately. 



3.3 Analytic reaction rate fits 
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3.3.1 Parametrization 



Reaction rates have been calculated for a temperature grid of 24 temperatures: 
T 9 =0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 
5.0, 6.0, 7.0, 8.0, 9.0, 10.0. For easy application in astrophysical investigations, 
all reaction types ((n/y), (n,p), (n,a), (p/y), (p,n), (p,a), (0,7), (a,n), (a,p), 
(7,n), (7,p), are fitted with the same parametrization 

exp (a + aiT^ 1 + a 2 Tg 1/3 + a 3 T g 1/3 + a 4 Tg 

+a 5 r 9 5/3 + a 6 lnT 9 ) , (16) 

with the seven open parameters ao — a$ and the stellar temperature Tg given 
in 10 9 K. This parametrization proves to be flexible enough to accommodate 
the different temperature dependencies of the various reaction types across the 
fitted temperature range of 0.01 < Tg < 10. Parametrizations of the present 
rates in the form used in Q and others can be obtained from the authors on 
request. 
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3.3.2 Parameters for the reverse rates 

The parameters for the reverse rates are not given explicitly but can easily be 
computed from the information in the tables. To calculate the reverse rate of 
the reaction i(j,o)m, i.e. the reaction m(o,j)i, Eq. ( |T^ ) is employed and the 
seven parameters a cv — ag cv for the reverse rate are determined as follows: 



rev 
"0 


— "0 1 


as tabulated 


a\ cv 


= ai - 


11.6045Q 


aT 


= a 2 




aT 


= a 3 




rev 
"4 


= a 4 




rev 
"5 


= a 5 




rev 
"6 


f a& 
\ a 6 


+ 1.5 (i) 
(ii) 



(17) 



The above relations are derived from Eqs. ( |l^ ) and (14), using Eq. jl^ ) and 
taking the logarithms on both sides. For the coefficient (i) applies 

when calculating a photodisintegration rate from a capture rate, case (ii) for all 
other rates. Finally, for the reverse reaction case, the value found by application 
of Eqs. fild[ ) and ( fT^ j has to be multiplied by the ratio of the partitio n func tions 
for residual and target nucleus Gi/G m . Examples are shown in Sec. 3.3.6] . 



3.3.3 Fit accuracy 

The flexibility of the fitting function makes it prone to numerical problems 
outside the calculated range at low temperatures, where the rates should be 
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close to zero. In some cases they tend to diverge strongly. This difficulty can be 
avoided by providing fit data at low temperatures additionally to the calculated 
values by appropriately extrapolating the rates to lower temperatures. This is 
achieved by either assuming s-wave capture for Tg < 0.1 for exoergic neutron 
capture reactions (Maxwellian averaged capture cross sections in the energy 
range 5 < E < 100 keV for targets along the line of stability can be found in 
another compilation JiTj]) or by considering proper Coulomb barrier penetration 
factors in the charged particle channels. Thus, both accuracy and flexibility can 
be ensured within a single parametrization. However, it has to be emphasized 
that the given parameterization is only valid within the temperature range of 
0.01 < Tg < 10., although many fits will show "proper" behavior down to lower 
temperature. Caution is advised when using derived (7,p) and (7,n) rates at 
the proton dripline (see below). For all cases, it is recommended to use the fits 
only down to the temperature given in the table. The temperatures of the 
validity of the fits are given in the tables for each reaction, to emphasize the 
importance of the given fit ranges. 

As a measure of the accuracy of a given fit, the quantity £ is quoted in the 
tables. It is defined by 

1 24 / f \ 2 

f-agG^) • < 18 > 

with r being the original rate value as calculated at each of the 24 temper- 
atures T 9 = 0.1, 0.15 . . . 10.0, and / the rate calculated from the fit at these 
temperatures. Contributions with r < 10~ 20 cm 3 s _1 mole -1 are neglected as 
lower accuracy at at low rates is inconsequential. Note that while a small value 
of £ is indicative of an accurate fit over the entire temperature range, large C 
generally signify deviations of the calculated from the fitted rate at the lowest 
temperatures only. 

The fit parameters are tabulated regardless of the validity of the statistical 



model of nuclear reactions in the given temperature range (see Sec. 2A). The 



estimated lower temperature limit of the validity of the statistical model, Tj^ 
is given separately for each rate in the tables. Below that limit the calculation 
of the rate by means of the statistical model may not be justified, although 
the fit to the calculated rate will still be accurate. At temperatures below 
the applicability limit, rates may be over-estimated and should be compared 
to calculations considering single resonance and direct reaction contributions. 
Especially close to the driplines, fits of reactions with low Q-value cannot be 
applied at low temperatures. Although the fit may be valid, it should not be 
used at low temperature because the statistical model will not be applicable 
anymore. 

3.3.4 Computed rate sets 

Two different sets of rates have been calculated. They differ in the mass model 
used, which enters into the computation of the separation energies and Q-values 
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as well as into the microscopic input to the level density calculation (see Sec, |2.3| ). 
One set was calculated employing the well-known FRDM mass model |4^|, which 
excellently reproduces masses and other ground state properties of nuclei close 
to stability. It is also the most comprehensive set across the nuclear chart. 

Recently, it was suggested that so-called shell quenching effects may arise 
for neutron-rich nuclei far off stability ]4£ , |50| . Fully microscopic calculations 
[ |5l[ |52| and experimental data JH^, [54], |55] , 56| indicate the weakening of nuclear 
shell gaps for neutron-rich nuclei. In the absence of microscopic calculations 
of the required nuclear properties for the whole nuclear chart, it is possible to 
phenomenologically include such quenching into existing mass formulae. This 
has been done in the ETFSI-Q model Q, based on the ETFSI-1 mass model @ 
|58| . However, it does not cover the full range of isotopes. Therefore, we provide 
the alternative sets of reaction rates obtained with the two mass models, so that 
the rates based on the FRDM can be used for large-scale studies and close to 
stability, and the rates based on ETFSI-Q for investigations concerning neutron- 
rich unstable isotopes and the r process. It has to be noted that one should 
refrain from mixing rates from the two sets as this will lead to inconsistencies 
and artificial effects in the results. 

3.3.5 Mass ranges of the tabulated fits 

Due to the extensive number of nuclear reactions in the considered mass range, 
we have to limit the printed version of our reaction rate fits. Full rate libraries 
both for reactions calculated with the FRDM and with the ETFSI-Q (as well as 
ETFSI-1) mass model can be obtained on-line or from the authors on request 
(see Sec. ||). In the printed version, only the FRDM set is given and no capture 
rates for reactions with negative Q value are shown. 

The full electronic versions of the tables available on-line include all reactions 
in the range 10 < Z < 83 (FRDM) and 24 < Z < 83 (ETFSI-Q). This amounts 
to 5369 (FRDM) and 4628 (ETFSI-Q) involved nuclei. The isotope ranges for 
which rate fits are available are given in Table A; for the FRDM, the mass range 
is also indicated by the heavy lines in Tables IA-IC. Rate fits are given for all n-, 
p-, and a-capture reactions and for those (n,p), (n,a), (p,n), (p,a), (a,n), and 
(a,p) reactions having positive Q-value. Reverse rates are not given expli citly 



but c an be computed by a two-step procedure as described in Sees. |3.3.2| and 
3.3. 6; . The stellar enhancement factors close to stability as well as the partition 



functions for all isotopes are given for a temperature grid of 24 temperatures: 
T 9 = 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 
4.5, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 . 

The printed Tables II— IV contain the following (calculated with the FRDM 
mass model): 

• (n,7) rate fits from stability to close to the neutron dripline in the range 
10 < Z < 49, and to N = Z + 59 for 50 < Z < 83 (Table II). 

• (n,p) and (n,a) rate fits around stability in the range 10 < Z < 83 (Table 
II). 
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(p/y)) (P) n ); (P) a )) ( a >7)> ( a ) n )) an d ( a )P) ra te fits from proton-rich nuclei 
to stability in the range 10 < Z < 50 (Tables III and IV). 

Reverse (endoergic) rates are not given explicitly but can be computed 
with the help of partition functions from the information given in the 
tables (see Sees. |3A2| , |3.3.6| ). 



Stellar enhancement factors /* at selected 16 temperatures are only quoted 
close to the valley of stability (Tables II— IV). 

• Partition functions for all involved isotopes are given at selected 20 tem- 
peratures in the range 0.1 < T 9 < 10 (Table V). 

An overview of the provided rates is given in Tables IA— IC, which show in 
detail for which n-, p-, and a- induced reactions rate fit parameters are available 
in Tables II— IV, respectively. 

3.3.6 Examples of use of tables 

This section is intended to help with interpreting the information given in the ta- 
bles. We give two examples for calculating the reaction rate for a given reaction 
and its inverse reaction at a temperature of Tg = 2.0. 

The first example is the capture reaction 35 Ar(p,7) 36 K. From Table III one 
finds a Q-value of Q — 1.666 MeV and the parameters ao =128.39, ai = —4.0033, 
a 2 = 137.67, a 3 = -276.87, a 4 =17.691, a 5 = -1.0728, a 6 =123.68. With the 
help of Eq. (|l6h one calculates Na (crv)* = 92.5 cm 3 s _1 mole _1 at Tg = 2.0. 
Because both T^w ani ^ ^k>w are considerably smaller than our temperature Tg, it 
is safe to assume that the statistical model is applicable and the fit to the rate is 
valid. In order to obtain the value for the reverse rate, one first has to determine 
the parameter values in the given parametrization. The parameter a cv = 151.83 
is given in the table. The remaining parameters are derived according to Eq. ( p"7| ) 
for a photodisintegration rate. This yields a^ cv = —23.3364, o,g cv = 125.18; all 
other parameters assume the same value as for the forward reaction. Using 
those, Eq. @ 

gives a value of A 7 ' =2.5xl0 8 s . This has to be multiplied by 
the ratio of the partition functions in order to obtain the valid rate factor for 
36 K( 7 ,p) 35 Ar: 



The values of the partition functions at Tg = 2.0 were taken from Table V. 
Note that for capture rates the procedure is always the same as described above 
regardless of whether it is an exoergic or an endoergic reaction. 

The second example we consider is the reaction 34 S(a,n) 37 Ar, again at 
T 9 = 2.0. It is not to be found in Table IV because of its negative Q-value. 
Therefore, one has to rely on Table II to calculate this reaction as the inverse 
reaction of 37 Ar(n,a) 34 S. The parameters found in Table II are ao =20.072, 
ai = -0.019613, a 2 = 1.8224, a 3 = -4.759, a 4 =0.56437, a 5 = -0.033893, 
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a e =1.7801, and Q = 4.63 MeV. Again, both T££ and are lower than 
the temperature of interest. The rate Na (<t«}* = 5.2 x 10 7 cm 3 s _1 mole _1 for 
37 Ar(n,a) 34 S results from the direct application of Eq. ([i"6|). In order to calcu- 
late the rate of 34 S(a,n) 37 Ar, the parameters are determined by application of 
Eq. @. This yields the value af v = -53.748448. The value of a r cv = 20.199 
is taken from the table and all other parameters remain the same as for the 
forward reaction. With Eq. ( fl6| ) one arrives at the rate value r' = 1.272 x 10~ 4 
cm 3 s _1 mole _1 at Tg = 2.0. Multiplying this by the appropriate ratio of parti- 
tion functions taken from Table V yields the final result 

N A (av)* = r'^r = r ,^0 = ^ ^ lQ _ 4 cm 3 g -i mole -i_ 
G345 1.0 



4 Summary 

Thermonuclear reaction rates for neutron-, proton- and a- induced reactions and 
their inverses have been calculated in the statistical model. All rates from the 
proton dripline to the neutron dripline for 10 < Z < 83 (Ne to Bi) have been fit- 
ted to a unique function with seven free parameters. Tables of these parameters 
are provided on-line for two sets of rates, calculated with input from two dif- 
ferent mass models. Furthermore, the stellar enhancement factors are given in 
order to facilitate comparison with experimental ground state rates. A printed 
subset of the on-line tables for the FRDM presented here shows fit parameters 
for (n/y) and (p,7) reactions from close to their respective driplines to stabil- 
ity, and for other n-, p-, and a- induced reactions with positive Q-values near 
stability. A prescription on deriving rates for inverse reactions with negative 
Q-values is given, as is a listing of the necessary partition functions. 

It should further be noted that only purely theoretical rates are given here 
which do not use any direct experimental information (except for nuclear masses 
and excited state information where available) . The methods to predict nuclear 
properties needed in the statistical model calculations are chosen to be as reliable 
as possible in order to retain predictive power. This is a compromise which 
may lead to locally enhanced inaccuracies but it emphasizes the importance of 
reliable predictions of rates far off stability. 

In real applications, these rates should be supplemented or replaced with 
experimental rates as they become available. Such a combination of theo- 
retical and experimental rates is provided, e.g., in the REACLIB compila- 
tion. Latest info rmation on the current vers ion of REACLIB can be found 
on the WWW at / http://ieJbl.gov/astro. htm\ /. Further details on the NON- 



S MOKER code and the cross section and reaction rate ca lculations are presented 
/ http: //quasar. physik.unibas.ch~tomrny/reaclib.htm\ /. Rates including fur- 



at 

ther mass models can also be obtained from the authors on request or directly 
at the latter URL. 



15 



Acknowledgements 

This work was supported in part by the Swiss National Science Foundation 
(grant 2000-053798.98) and the Austrian Academy of Sciences (APART). T. 
R. is a PROFIL fellow of the Swiss National Science Foundation (grant 2124- 
055832.98). We thank S.E. Woosley and R.D. Hoffman for discussions. We also 
want to thank the Consulting Editor Wilfried Scholz for helpful comments and 
help in the preparation of the final manuscript. 

References 

[1] M. Arnould, Astron. Astrophys. 19, 92 (1972). 

[2] J. Holmes, S. Woosley, W. Fowler, and B. Zimmerman, Atomic Data Nucl. 
Data Tabl. 18, 305 (1976). 

[3] S. Woosley, W. Fowler, J. Holmes, and B. Zimmerman, Atomic Data Nucl. 
Data Tabl. 22, 371 (1978). 

[4] D. G. Sargood, Phys. Rep. 93, 61 (1982). 

[5] F.-K. Thielemann, M. Arnould, and J. W. Truran, in Advances in Nuclear 
Astrophysics, edited by E. Vangioni-Flam et al. (Editions frontiere, Gif sur 
Yvette, 1987), p. 525. 

[6] J. J. Cowan, F.-K. Thielemann, and J. W. Truran, Phys. Rep. 208, 267 
(1991). 

[7] L. Wolfenstein, Phys. Rev. 82, 690 (1951). 

[8] W. Hauser and H. Feshbach, Phys. Rev. 87, 366 (1952). 

[9] T. Rauscher, F.-K. Thielemann, and K.-L. Kratz, Phys. Rev. C 56, 1613 
(1997). 

[10] C. Mahaux and H. Weidenmiiller, Ann. Rev. Part. Nucl. Sci. 29, 1 (1979). 

[11] T. Rauscher and F.-K. Thielemann, in Stellar Evolution, Stellar Explosions 
and Galactic Chemical Evolution, edited by A. Mezzacappa (IOP, Bristol, 
1998), p. 519. 

[12] S. Goriely, in Nuclei in the Cosmos V, edited by N. Prantzos and S. Haris- 
sopulos (Editions Frontieres, Paris, 1999), p. 314. 

[13] R. B. Firestone, in Table of Isotopes, 8th ed., edited by V. S. Shirley (Wiley, 
New York, 1996). 

[14] P. Moller, J. R. Nix, and K.-L, Kratz, Atomic Data Nucl. Data Tabl. 66, 
131 (1997). 



16 



[15] J. Jeukenne, A. Lejeune, and C. Mahaux, Phys. Rev. C 16, 80 (1977). 

[16] S. Fantoni, B. Friman, and V. Pandharipande, Phys. Rev. Lett. 48, 1089 
(1981). 

[17] C. Mahaux, Phys. Rev. C 82, 1848 (1982). 

[18] L. McFadden and G. R. Satchler, Nucl. Phys. 84, 177 (1966). 

[19] P. Mohr et al, Phys. Rev. C 55, 1523 (1997). 

[20] E. Somorjai et al, Astron. Astrophys. 333, 1112 (1998). 

[21] G. R. Satchler and W. G. Love, Phys. Rep. 55, 183 (1979). 

[22] T. Rauscher, in Proc. XXVI Int. Workshop on Gross Properties of Nu- 
clei and Nuclear Excitations, edited by M. Buballa, W. Norenberg, J. 
Wambach, and A. Wirzba (GSI, Darmstadt, 1998), p. 288. 

[23] T. Rauscher, in Nuclear Astrophysics, Proc. 9th Ringberg Workshop, edited 
by W. Hillebrandt and E. Miiller (MPA, Garching, 1998), Vol. MPA/P10, 
p. 84. 

[24] J. M. Blatt and V. Weisskopf, Theoretical Nuclear Physics (Wiley, New 
York, 1952). 

[25] E. Lipparini and S. Stringari, Phys. Rep. 175, 103 (1989). 
[26] W. D. Myers et al, Phys. Rev. C 15, 2032 (1977). 
[27] M. Danos, Nucl. Phys. 5, 23 (1958). 

[28] F.-K. Thielemann and M. Arnould, in Proc. Int. Conf. on Nucl. Data for 
Sci. and Techn., edited by K. Bockhoff (Reidel, Dordrecht, 1983), p. 762. 

[29] C. M. McCullagh, M. L. Stelts, and R. E. Chrien, Phys. Rev. C 23, 1394 
(1981). 

[30] S. Raman, O. Shahal, and G. G. Slaughter, Nature 331, 505 (1981). 

[31] S. G. Kadmenski, V. P. Markushev, and V. I. Furman, Sov. J. Nucl. Phys. 
37, 165 (1983). 

[32] J. Kopecky and M. Uhl, Phys. Rev. C 41, 1941 (1990). 

[33] T. Rauscher, F-K. Thielemann, J. Gorres, and M. C. Wiescher, Phys. Rev. 
C, submitted (2000). 

[34] J. J. M. Verbaatschot, H. A. Weidenmiiller, and M. R. Zirnbauer, Phys. 
Rep. 129, 367 (1984). 

[35] J. Tepel, H. Hoffmann, and H. Weidenmiiller, Phys. Lett. 49B, 1 (1974). 



17 



[36] E. Gadioli and P. E. Hodgson, Pre-Equilibrium Nuclear Reactions (Claren- 
don Press, Oxford, 1992). 

[37] S. Ezhov and V. Plujko, Z. Phys. A 346, 275 (1993). 

[38] J. Thomas, M. Zirnbauer, and K.-H. Langanke, Phys. Rev. C 33, 2197 
(1986). 

[39] A. Gilbert and A. Cameron, Can. J. Phys. 43, 1446 (1965). 

[40] D. J. Dean et al, Phys. Rev. Lett. 74, 2909 (1995). 

[41] V. Paar and R. Pezer, Phys. Rev. C 55, R1637 (1997). 

[42] P. Moller, J. R. Nix, W. D. Myers, and W. J. Swiatecki, Atomic Data Nucl. 
Data Tabl. 59, 185 (1995). 

[43] J. M. Pearson, R. C. Nayak, and S. Goriely, Phys. Lett. B 387, 455 (1996). 

[44] G. Audi and A. H. Wapstra, Nucl. Phys. A595, 409 (1995). 

[45] T. Rauscher et al, Phys. Rev. C 57, 2031 (1998). 

[46] W. A. Fowler, Quarterly Journ. Royal Astron. Soc. 15, 82 (1974). 

[47] Z. Y. Bao et al, Atomic Data Nucl. Data Tabl., in press (2000). 

[48] W. A. Fowler, G. E. Caughlan, and B. A. Zimmermann, Ann. Rev. Astron. 
Astrophys. 5, 525 (1967). 

[49] K.-L. Kratz et al, Ap. J. 403, 216 (1993). 

[50] F.-K. Thielemann et al, Nucl. Phys. A570, 329c (1994). 

[51] B. Chen et al, Phys. Lett. B 355, 37 (1995). 

[52] J. Dobaczewski et al, Phys. Rev. C 53, 2809 (1996). 

[53] T. R. Werner et al, Phys. Lett. B 333, 303 (1994). 

[54] P. Hoff et al, Phys. Rev. Lett. 77, 1020 (1996). 

[55] T. Kautzsch et al, Phys. Rev. C 54, 2811 (1996). 

[56] C. T. Zhang et al, Z. Phys. A 358, 9 (1997). 

[57] Y. Aboussir, J. M. Pearson, A. K. Dutta, and F. Tondeur, Nucl. Phys. 
A549, 155 (1992). 

[58] Y. Aboussir, J. M. Pearson, A. K. Dutta, and F. Tondeur, Atomic Data 
Nucl. Data Tabl. 61, 127 (1995). 



18 



5 Explanation of Tables 



Table IA: Neutron-Induced Reaction Rates Available in Ta- 
ble II 

This is an overview of which neutron-induced reaction rates are available in the 
printed and the online versions. The full lines delimit the range of rates in the 
electronic version as given in Table A for the FRDM. The entries at a single 
neutron and proton number specify the reactions on the given target nucleus 
listed in the printed Table II. Only reactions with positive Q- va lue are shown. 
In addition to the marked rates, their reverse rates (with negative Q- value) can 
be inferred from the information in Table II as explained in Sec. |3.3,2| . The 
reactions are denoted as follows: 

G (n, 7 ) 

P (n,p) 

A (n,a) 

The box at the lower left corner gives the location in the Z, N plane of the 
final nucleus relative to the target nucleus for (n,7), (n,p), and (n,a), thereby 
specifying also the inverse reaction fits derivable from Table II. 

Table IB: Proton-Induced Reaction Rates Available in Ta- 
ble III 

Same as Table IA but for proton-induced reactions. The marked reactions 
correspond to the entries in Table III. The reactions are denoted as follows: 

G (p,7) 
N (p,n) 
A (p,a) 

The box at the lower left corner gives the location in the Z, N plane of the 
final nucleus relative to the target nucleus for (p,7), (p,n), and (p,a), thereby 
specifying also the inverse reaction fits derivable from Table III. 

Table IC: Alpha Particle-Induced Reaction Rates Available 
in Table IV 

Same as Table I A but for a-particle induced reactions. The marked reactions 
correspond to the entries in Table IV. The reactions are denoted as follows: 

G (a l7 ) 

N (a,n) 
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P (a,p) 



The box at the lower left corner gives the location in the Z , N plane of the 
final nucleus relative to the target nucleus for (a, 7), (a,n), and (a,n), thereby 
specifying also the inverse reaction fits derivable from Table IV. 



Table II: Neutron-Induced Reaction Rates 

Fits to stellar rates Na (cu)* for (11,7), (n,p), (n,a) reactions, calculated includ- 
ing masses from the FRDM. The rates in cm 3 mole -1 s _1 are computed by 
the use of Eq. (|ltf), with the temperature given in units of 10 9 K. The fits are 
valid in the temperature range Tj^ < Tg < 10, with given in the table. It 
should be noted that while the fit may still be formally valid and accurate, the 
application of the statistical model may not be justified at low temperatures. 
An estimate for the applicability of the statistical model is given by Xj^. The 
following information is provided: 

Target Reaction target 

Reaction Reaction type and final nucleus 

Q Reaction Q- value 

Ji Target ground state spin (same as in Table V) 

Jf Final nucleus ground state spin (same as in Table V) 

Tj™ Estimate of the lower temperature limit for the applicability of the 

(Hauser-Feshbach) statistical model; "n.c." indicates that the limit 
was not calculated for the given reaction. 

Lower temperature limit for the fit; usually 0.01. Note that the fits 
give extrapolated rates below Tg = 0.1, which may be less accurate, 
especially if they are very small. 

Dev Fit accuracy £ (Eq. (|l8|)) 

SEF If the field is blank, the seven fit parameters below are followed 

by the stellar enhancement factors (SEF) /* (Eq. (|ll])) at the 16 
temperatures given in the head of the table. 

A value of 1 indicates that all SEF are unity; no SEF are printed. 
A value of indicates that no SEF were calculated. 

do. . . ag Seven fit parameters for the forward rate 



a T Q V First parameter for the reverse rate fit (see Sec. 3.3.2 

J.SEF Temperatures at which the SEF were calculated 
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Table III: Proton-Induced Reaction Rates 

Same as Table I for the reaction types (p/y), (p,n), (p,a). 

Table IV: Alpha Particle-Induced Reaction Rates 

Same as Table I for the reaction types (0,7), (a,n), (a,p). 

Table V: Partition functions 

Partition functions of isotopes for various temperatures calculated with a level 
density making use of FRDM input. Included are only those partition functions 
for nuclei involved in the reactions given in Tables II— IV. 

Nuc Isotope for which the partition functions are tabulated. 

Tq Temperature (in 10 9 K) at which the partition functions have been 

calculated. 

P A value of 1 indicates that all partition function are unity; no par- 

tition functions are then printed explicitly. 

Spin Ground state spin of nucleus, either from experiment |Q or from 

theory JH). 

Partition Functions: Partition functions normalized to the ground state (Eq. (p^)) 
for the 20 temperatures specified in the table header. 
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Table 1: This table lists the isotope range of the full rate tables which are 
available electronically. Given are the charge number Z of a target and the 
lower and upper limits iV m i n and N max of the neutron number in the isotopic 



chain. 





FRDM 


ETFSI-Q 




FRDM 


ETFSI-Q 


z 


AT min 


Nmax 




Z 


JVmin 


N 

mm 




Af max 


8 


5 


10 




47 


41 


113 


41 


111 


g 


5 


28 




48 


42 


115 


42 


112 


10 


5 


31 




49 


43 


117 


43 


113 


H 


6 


33 




50 


44 


119 


44 


114 


12 


7 


35 




51 


46 


121 


46 


115 


13 


8 


38 




52 


47 


124 


47 


124 


14 


8 


40 




53 


48 


126 


48 


126 


15 


8 


42 




54 


49 


128 


49 


128 


16 


8 


44 




55 


51 


130 


51 


130 


17 


9 


46 




56 


52 


133 


52 


132 


18 


9 


49 




57 


53 


135 


53 


133 


19 


10 


51 




58 


55 


137 


55 


134 


20 


10 


53 




59 


56 


139 


56 


135 


21 


11 


55 




60 


58 


141 


58 


136 


22 


12 


58 




61 


59 


144 


59 


137 


23 


13 


60 




62 


61 


146 


61 


138 


24 


14 


62 


18 62 


63 


62 


148 


62 


139 


25 


15 


64 


18 64 


64 


64 


150 


64 


150 


26 


16 


66 


19 66 


65 


65 


153 


65 


152 


27 


17 


69 


19 67 


66 


67 


155 


67 


154 


28 


18 


71 


20 68 


67 


69 


157 


69 


155 


29 


19 


73 


21 69 


68 


70 


159 


70 


156 


30 


21 


75 


22 70 


69 


72 


161 


72 


157 


31 


22 


77 


23 71 


70 


73 


164 


73 


158 


32 


23 


80 


24 72 


71 


75 


166 


75 


159 


33 


24 


82 


25 73 


72 


77 


168 


77 


160 


34 


25 


84 


26 84 


73 


78 


170 


78 


161 


35 


26 


86 


27 86 


74 


80 


173 


80 


162 


36 


27 


88 


28 88 


75 


81 


175 


81 


163 


37 


29 


91 


29 89 


76 


83 


177 


83 


177 


38 


30 


93 


30 90 


77 


85 


179 


85 


179 


39 


31 


95 


31 91 


78 


87 


182 


87 


182 


40 


32 


97 


32 97 


79 


88 


184 


88 


184 


41 


33 


99 


33 99 


80 


90 


186 


90 


186 


42 


35 


102 


35 102 


81 


92 


188 


92 


188 


43 


36 


104 


36 104 


82 


93 


191 


93 


191 


44 


37 


106 


37 106 


83 


95 


193 


95 


193 


45 


38 


108 


38 108 


84 


98 


193 


98 


193 


46 


40 


110 


40 110 


85 


101 


195 


101 


195 
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The remaining tables can be found on the ADNDT 
server 

or at http : / / quasar . physik . unibas . ch/~tommy/ adndt . html . 
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